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Correction terms generated in the correlator analysis due to multiplicity-dependent 

66 ■ 

, observable mean are investigated. A procedure for subtraction of such terms from 

o ■ 

^ ' calculated correlator estimates is suggested and the obtained results are discussed. 

I. INTRODUCTION. 

Behavior of strongly interacting matter under extreme conditions attracts interest of 
O ' physicists for a considerably long time. Motivated partially also by attempts to explain 

^ ! n 

^ . interesting cosmic-ray observations [1] it has been assumed, that in relativistic collisions 



O 



O 



of hadrons or nuclei a short-living system containing a super-dense thermally equilibrated 
matter can be created in laboratory. Intriguing results of experimental groups at CERN and 



in 

> 

00 

^ i BNL show us that such an effort is a challenging task and interpretations of the obtained 



experimental information are still vividly discussed. 

It has been suggested [2I] that calculation of the correlator ratios may allow us to estimate 
quantitatively a degree of thermalization of high-multiplicity thermodynamical systems cre- 
ated in relativistic collisions of hadrons or nuclei. This suggestion has been met with a 



. considerable interest 



and recently a procedure for the evaluation of higher-order correla- 



tors in a reasonable computer time has been presented {4]. 

On experimental side correlator Cf* has been evaluated using transverse momenta of 
charged particles created in relativistic Au+Au collisions at RHIC [s^ and centrality depen- 
dence of such correlator has been studied. However, it seems that calculation of correlators 
from real experimental data requires a more detailed approach. For example, if the observ- 
able mean of quantity (e.g. (pt)) used for calculation of correlators depends on measured 
multiplicity, correlators can be systematically shifted and obtained results biased. 
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In the following sections we study quantitatively the influence of multiplicity-dependent 
observable mean on correlators and describe a simple procedure allowing one to subtract 
the introduced systematical bias from the calculated correlator estimates. 



II. A SIMPLE CORRECTION TO C2 CORRELATOR 

Let us start with global 2-particle correlator C2^'^^ calculated from measured character- 
istics Xi (e.g. transverse momentum pt) of observed charged particles with multiplicity n. 
For the sample A containing Ngy events one has (see Eq.(21) in and also Eq.(2) in 5]) 

Mn,x) _ 1 v^i^^ {Xj - {xG))iXj - {xg)) f.. 

where Nij = n\/{n — 2)\ is the number of all particle pairs in a given event with indexes i j 
(symbol [Nij] reminds double-sum) and (xg) is the global mean of particle characteristics 
Xi (e.g. transverse momentum pt) over all events and particles used for the analysis. 

Because correlators are not sensitive to number of particles n selected for their calcu- 



lation (see Section 3 and Tab. 3 in [4]) one can select a fixed number of tracks n < for the 
correlator evaluation in the sample of events with multiplicities ni < Uk < n2- 

It can be easily shown that if observable mean x(nfc) in events with observed multiplicity 
Uk depends on event multiplicity Uk then calculated value of correlator C2 for sample A is 
shifted by the amount A2 which depends on the width of multiplicity interval (ni,n2) and 
on properties of multiplicity-dependent observable mean x{nk). 

if ^ ; C^""'" = C*™' + A2 (2) 

In order to demonstrate this let us divide event sample A into subsamples Al and A2 
with multiplicities (ni,?T,ia) and {nia,n2) and let us assume (see Fig.l) that observable 
mean x{nk) increases with multiplicity in step-like way: x{ni < < riia) = xai and 
x{nia < "".fc < ^2) = xa2- For correlator Cg"'^'* calculated for event sample A one can write 

[^-■1 (x. - {xg)){x, - {xg)) , ^ (a;,, - {xG)){xi - {xg)) 



^A=A1+A2 _ ^ 



^ Iv'. ^ TV- ■ 

Nai i^j ^^'J Na2 iy^j 



(3) 



It is clear that in subsample Al with observable mean xai the correlator C2 is calculated 
using an overestimated global mean {xg) which introduces a static shift of global mean 
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Axai = xai — {xq) in the correlator calculation. Since (xj — {xg)) = {xi — xai + Axyn) one 
obtains for events within sub-sample Al (and similarly A2): 

{xi - {xG)){xj - {xg)) = {xi - XAi){xj - Xai) + AxAi[xi - Xai + Xj - xai] + (As^i)^ (4) 

A substitution into Eq.(3) using J^NAk Sr=i(^« ~ ^Ak) = (for k=l,2) gives 



^Al+A2 ^ Nai ■ Cj' + Na2 ■ Cf ^ Nai ■ {AxAl? ^ Na2 ■ (Axa2 



where the correlator for sub-sample Al (and similarly for A2) is 



Nai 



E 



] {xi - XAi){xj - Xai) 



Nai 



N. 



(5) 



(6) 




Fig.l: Calculation of C2 correlator using sub-samples Al and A2. 

Assuming Nai = = Nev/2 and |Aa;A2| = |Axyii| (see Fig.l) one obtains a simple 
relation 



2 4 

where A{xg) = xa2 — xai = 2Axa2- If dynamics of particle production in sub-samples 
Al and A2 is identical then ~ and one can express the exact "true" correlator 
C^^^^ = {C^^ + Cf)/2 for the event sample A as 



a 



true 



a 



calc (A(xg))' 



4 



(8) 



Correction term (A(xg'))^/4 can be calculated and subtracted from the correlator C2°'^'^ 
which is evaluated using global observable mean {xg) of the whole sample A = Al + A2. 
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Increasing the number of sub-samples Al, A2, . . . , gives (in general) a more precise cor- 
relator estimate. Using 4 subsamples of the event sample A — Al + A2 + A3 + AA one 
obtains 

' Nai + Na2 + Nas + Na4 4 ^ ^ ^ ^ ' 

if multiplicity dependence of observable mean x{nk) increases in 4 steps: x{nk&Ai) = xai, 
x{nkeA2) = XA2, x{nkeA3) = xas, xiukeA^) = xa4 (see Fig.2). 




Fig.2: Calculation of C2 correlator using subsamples Al + A2 + A3 + A4. 



The goal of the next section is to show how to evaluate correction terms for any given 
dependence of observable mean xaink) on observed event multiplicity n^. 



III. CORRECTIONS TO C2 IN GENERAL CASE 

Let us consider a general multiplicity-dependence of observable mean x{nk). After sub- 
stitution (xg) — ^ (xg) — x{nk) + x{nk) in Eq.(l) correlator C2 is 

(jA ^ ^calc i^i - ^(^fc) + ^XG{nk)][Xj - xjuk) + Axgink)] ^^^^ 

where AxGijik) = x{nk) — (xg)- Here observable mean x{nk) of events with measured 
multiplicity Uk is defined as 



X[ 



N=l 



^^)=E E?f/^n, (11) 



=1 ^ 



where N^^ is total number of events in sample A with measured multiplicity Uk and n 
denotes number of particles (in each event) used for the calculation of correlators. Since 
Z]jv=iSr=i(^i ~ x{nk)) — the expression becomes 

^caic_ 1 ^ (xj - x{nk)){xj - xjuk)) , ^ {AxG{nk)f , . 

^2 -]^2.1^ + ^ — N ^ ^ 
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The first term on the right-hand side of Eq. (fT2|) is a correlator calculated using a 
multiplicity-adjusted observable mean x{nk). We call this "true" correlator Cg*'''^'^ which 
coincides formally with global correlator defined by Eq.([T]) evaluated for events with fixed 
multiplicity Uk- Such correlator does not contain correction terms of type given by Eq.(7) 
due to zero width {n2 = rii = Uk) of multiplicity interval of events being analyzed. True 
correlator 

^trne ^ J_ ^ (x— x(nfc)) (Xj - x(nfc)) ^^^^ 

^ev i-^j Nij 

is also a limiting case of the first term on the right-hand side of Eq.([9]): 

2 k->An Y!l=lNi ^ ' 

where An = n^ — is the maximal number of multiplicity sub-intervals in the sample A. 
Using a more explicit notation one has 

Cr\x{nk)) = Cr'^iixa)) - E ^^^^^^ ~ ^^^^^^ (15) 

/V * en 

ev 

where correction term J2{x{nk) — (xg))^ /^^ev can be further expressed as 

Here is the number of events with multiplicity in sample A and probability density 
P{nk) for events with multiplicity Uk is P{nk) = N^^/ J2 ^^n^- Based on Eq. ffTSj) one can 
find "true" correlator for the event sample A by subtracting the correction term defined 
by Eq.(fT6l) from global correlator Cg'^^'"^ = C^((xg)) (calculated using the global observable 
mean (xg) = J2xG{i^k)Nn^/Nev)- To summarize our result in analytical form we express 
"true" 2-particle correlator as 



^ true 



'2 



x{nk)) = Cr'^((a;G)) - / P{nk) ■ [x{nk) - {xG)]'dnk (17) 



where meaning of P{nk), (xg) and x{nk) has been described in the text above. Assuming 
linear multiplicity dependence of observable mean {x{nk) = Xq + k ■ Uk) and constant P{nk) 
distribution one evaluates integral in Eq. (fT7|) as k'^{n2 — ^1)^/12 = (A(xg))^/3. This is in 
agreement with Eq.Q suggesting corrections [1 + 1/4 + 1/16 +■ ■ ■ ](A(2;c;))^/4 if number of 
subsamples {A^} is iteratively doubled. (Note that T,n=o 1/2^" = 4/3.) 



6 



IV. CORRECTIONS FOR Ck CORRELATORS 



Results obtained for C2 correlator can be generalized for higher order correlators. Let us 
define i^-th order global correlator Ck 



1 



(18) 

where Ni^i2...ij^ = n\/{n — K)\ is the number of particle i^-plets {«i,^2---^x} niade from 
n < rik particles in each event (n^ is total multiplicity of a given event). After substitution 
{xg) =x{nk) — Axcini:) (where AxG'(?7,fc) =x{nk) — {xg)) a compact expression can be found 



^caic..- .^ (AxG(nfc ))^-^ Cr^(xG(n,)) f 

A = iVe« 



where C^^^^^lxlni^)) are A-th order "true" correlators (defining Cg*™'^(a;(nfc)) = 1) formally 
obtained by replacing (xg) in Eq. (fT8|) by x{nk). In particular, for = 3 one has 

Nev Nev 

since C^'^^'^{x{nk)) = due to definition of x{nk). For constant probability P{nk) and 
for function AxG{nk) antisymmetric around average (ra^) term {AxG{nk)Y becomes zero. 
Moreover, if C2^'^'^{x{nkj) is a constant function of Uk one obtains C^^^^{{xg)) = C*3™°(a;(nfc)). 
Under the same assumptions for N = 4 one has 



and for = 5 



(AxG(nfc))^ ^ g (AxG(n,))2 C2*'^-(x(nfe)) 



A^. 



Cr'^(x(n,)) = Cr'=((xG)) + 10E 



(AxG(nfc))2Cf"^(x(nfc)) 



A",, 



(21) 



(22) 



where C;[^'^"'^(a;(nfc)) are the A-th order true correlators defined as 

1 



Aij_i2---«A 



(23) 



We have thus obtained expressions for differences between C^^'^ = C^^^{{xg)) correlators 
calculated using global observable mean {xg) and true correlators C^^'^ = C^'^'^{x{nk)) 
evaluated using multiplicity-adjusted observable mean xijik): 



C 



K 



C^^"^+ A 



K 



(24) 
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Correction term A2 according to Eq. (fT7|l is 

A, = E^^^#^= r>W[^H-(^G)Frfn (25) 

AT ^^ev -^rii 

and from Eq. fl2Tl) one has 

A4 = / Pin) [x{n) - {xg) fdn + 6 / P{n) C2*™"(x(n)) [x{n) - {xg) fdn (26) 
For P{nk) symmetrical around (n^) term A3 —>■ 0. However, correction A5 remains non-zero 

A5 = 10 / Pin) "^(x(n)) [x(n) - (xg) ? dn (27) 

J ni 

For probability distribution V{nk) asymmetrical around (n^) one obtains from Eq. (l20|) 

rn2 rn2 

A3 = / V{n) [x{n) - {xg) fdn + / P(n) C2*™"(x(n)) [x{n) - {xg) ]dn (28) 

Jni Jni 

In a real calculation one does not have to calculate explicitly the multiplicity-dependent 
observable mean x{nk) for each multiplicity n^. It is enough to split event sample A into 
reasonable number of subsamples Ai, A2, ■ ■ ■ and calculate mean x[yl7v] for subsamples 
An. Fitting x[ylAr] values with smooth function X{nk) gives approximation for x{nk). 

One can evaluate "true" correlators also directly using multiplicity-adjusted observable 
mean x{nk) for subsets of events with fixed multiplicities n^. This method has been utilized 
in pt correlations analysis [5]. However, there is still another systematical effect which 
should be accounted for if one wants to obtain correlator values free from contributions due 
to multiplicity-dependent observable mean x{nk). We discuss this issue in the next section. 

V. FIXED-MULTIPLICITY CORRECTIONS TO CORRELATORS 

Let us consider now the situation when multiplicity of particles nt in a given event is not 
known precisely and instead of n*°* (total number of particles produced in a given event) a 
multiplicity of tracks hk ~ 72*°*/^ is measured in a detector. In events with a given fixed 
measured multiplicity there will be fluctuations of corresponding n^°* around the average 
value (n*°*) ~ nk-^ (where ^ > 1 is a real number corresponding to the detector acceptance). 
If observable mean of quantity under study does not depend on multiplicity {x{n) = const) 
there is no influence on correlator values from these fluctuations. 



8 



However, if observable mean x{n) depends on multiplicity additional corrections A„ to 
calculated correlator values appear due to n*°* fluctuations at fixed measured hk (which 
generate fluctuations of observable mean x{n\°^) values at given hk). 

Contributions A„ are the fixed-multiplicity corrections to correlators. They can influence 
results and interpretations of the correlator analysis if they are not accounted for. 

For large enough measured multiplicities fluctuations of measured multiplicity hk at 
given fixed total multiplicity n^°^ are close to Gaussian (see Appendix) with probability 
distribution P{hk\n^°^): 

P{nk\nn = . , (29) 



where {h^) ~ ^"^^ '^fik ~ ^' \/(^*°V'C) (see Appendix). One can express probability 

of n*°* fluctuations P(n*°*|nfc) at given fixed measured using Bayes' theorem [gI. For 
constant P(n*°*) = const probability distribution one obtains (see Appendix) 

P(n*°*|nfc) = ^= (30) 



tot 



where cxtot = ^ ■ <^hk- Assuming linear approximation x(n*°*) = xq + k ■ n^°^ for multiplicity- 
dependent observable mean x(n*°*) in Scrtot vicinity of ^ hk ■ ^ one obtains fixed- 
multiplicity correction term A2 from Eq. fl25|l as 



A, = !p{n'°'\hk) [x{n'^') - x{{n'°%,)fdn'°'=~e f ^''''> " -dn (31) 

This gives a simple result 

A, = ~e-aUO (32) 

where o"tot(0 = '^fikC = CV^k {d') (see Appendix). For slope parameter k obtained from 
approximation x{hk) = Xq + k ■ h^. (measured experimentally) one has k = k/C, and thus 
A2 = k'^hk{(y)'^ where {a) 1 is to be obtained from MC simulation. Using similar 
arguments correction term A4 can be calculated from Eq. (!26|) as 

A4 = 3 fc^ ■ at, (0 + 6 P ■ al, (0 C^''^'' (33) 

Acceptance parameter ^ disappears {k = k/^ and (Xtot = f^nfeO- ^4 = ^ fc^o"!^ + Qk'^a\^C2^^ 
(where = y/uk (a) see Appendix). Fixed-multiplicity correction terms thus depend only 
on experimentally measurable quantities. 
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One might be tempted to imply A3 — > based on symmetrical probability distribu- 
tion given by Eq.( !30l) . However, fluctuations of (n*°* — values can be significantly 
asymmetrical for rik small enough (see Fig.3 for ?t,*°* fluctuations at fik = 10). 




Fig.3: Fluctuations of nf* for hk=W, ^=8.51 and nfc = 112, ^=2.44 (MC simulation). 

In this case correction A3 can be evaluated using Eq. fl28p . Symmetrical P(n*°*|rifc) given 
by Eq.(l29]) yields A3 = and from Eq.(l27D one obtains A5 = 10 k'^al^C^''''' . 

VI. MC SIMULATIONS 

In order to verify behavior of crtot(0 ^"^^ (^n^i'^k) a simple MC simulation has been 
performed: Events with total multiplicities n^°* G (50, 2000) have been generated with 
constant probability P{n^°^) = const. In each event, rapidities i/i were assigned to 
particles according to the bell-shaped rapidity distribution and number of observed particles 
hk found in the selected acceptance range {—ya, Va) has been determined. 

Two-dimensional histogram H2(nfc, n^°*) filled with pairs of obtained numbers n)^^ and fik 
is shown in Fig. 4. Projection histograms H„tot(nfc) and Hf,^(n^°*) which are proportional to 
probabilities P{nk\n^°^) from Eq. fHHj) and P(n*°*|nfc) given by Eq. flSlI) are shown for ^ ^ 4.26. 

Widths <7tot{C,, nl) and crfj^(n*°* = ^ ■ nl) have been obtained from Gaussian fits of the 
projection histograms. Width cTfij. has been found to follow c\frik dependence with c ~ 0.9 
in agreement with Eq. (H9l) . Ratio of widths o"tot(^, ""-fc) and o"nj.(r;,*°*) at ra*"*^ 'C'^fc appeared 
to be constant and equal to o"tot(^, '"■*°*/0/'^»ifc(^l°*) ~ ^ as expected from Eq.( l5Tl) . 

Simulations for P(n*°*) 7^ const have also been done and shift 6n^°^ = (^*°*)nj. — fik'C 
(described in the next secion) has been observed. 
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Fig.4: Histograms Rn^{n]^^), H„t^ot(nfc), H2(nfc,n^°*), fitted an^^ = c- .Jn^ and ratio (Xtot/o-fifc- 

VII. FIXED-MULTIPLICITY CORRECTIONS IN REAL EXPERIMENTS 

In a general case probability P{n^"^) ^ const (usually P(n*°*) decreases with The 
situation is slightly more complex in this case and there appears a shift of average {n^°^)nk 
(evaluated for events with fixed measured n^) relative to (fi*°*) = n^-^ expected in the case 
of P(n*°*) = const. We will show that the shift 

5n'^' = {n'"%^ -fik-i (34) 

does not influence fixed-multiplicity correction terms for correlators significantly. Let us 
assume P(n*°*) decreasing with n*°* as ao — ai ■ n*°* which gives (denoting axja^ = ui) 
normalized prior probability P(n*°*) in the form^ 

P(n*°*) = 2cj(l-cj-n*°*) (35) 

For a given fixed n*°* the average (n^) = n*°*/^ and if (n^) is large enough fluctuations 
of fik are close to Gaussian according to Eq.l lSTI) . However, average value (n*°*) for a given 
fixed measured hk is 

{n'°')r^^ = Jn''''P{n''''\fik)dn'°' (36) 

where P(n*°*|nfc) is a conditional probability distribution of having event with total multi- 
plicity ?T,*°* (for events with fixed measured multiplicity n^). Using prior probability P(n*°*) 
given by Eg. (135!) P(n*°*|nfc) can be found analytically if + 3 an^C, < l/Co (see Appendix) 



^ This probability distribution implies rt*°* < 1/uj. 
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e 



-in 



tot 



' ^ "fc (1 — n^^ct;) 



P(n'°'\nk) = ^ ^— ^ (37) 

and Eg. (15^ yields 

=i-n,- (38) 

One can also define effective observed multiplicity {nk)nk = (^*°*)nfe/^ for events with 
fixed measured hk 

(^fc)nfe = Uk- .^^,'^7'- = ^k- ^ ^^I"' (39) 

1 - UkiuJ 1 - rikuj 

Using slope parameter uj determined from P{nk) = 2uj{l —uo-hk) distribution (accessible 
experimentally) one has uj = u/C, comparing to Eg. (1351) and parameter ^ disappears from 
Eg. (139!) . Displacement 6hk = {nk)n^ — hk of the effective multiphcity {nk)hk relative to 
fixed measured multiplicity hk can be relevant e.g. in elliptic flow analysis where effective 
participant eccentricity e(^n^) is compared to elliptic flow strength V2{hk). 

Fixed-multiplicity correction term A2 for prior probability P(ra*°*) given by Eg. flSSl) can 
be evaluated analogously to Eg. (1311) 

e 

which for shifted {n^°^)n^ given by Eg. (1381) and P{n^°^\hk) given by Eg. (1371) yields 



^2=7?/ P(n"\At)ln"" - {n")ftJ'A" (40) 



Analytical expression for A4 can be found using Eg.( l26i) . Assuming (72™*^= one has 



At = 3 k'al 



1 - 



3A:MJ1-2.;VIJ (42) 



If probability P(n*°*) = const correction term A3 — > for hk large. For 7^ which 
means P(n*°*) 7^ const correction A3 can be found using Eg. fl28|) : 

2k^u^al , , . 

[l-nkUjY 

For sample of events with constant P(n*°*) distribution (which means constant measured 
P{hk) distribution) one has a; ^ 0; {n^°^)n^ =n^. ^- ifik)^^ = hk and A^ — > A„. 

Since acceptance parameter ^ is not present in Eg. fl41|42ll43p fixed-multiplicity correction 
terms can be determined from experimentally accessible guantities: k"^, uY and a\ . 
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VIII. RELATIONS BETWEEN CORRELATORS 

Let us assume that mean transverse momentum pt of particles in events with fixed total 
multiplicity n*°* (e.g. selected from the output of a MC event generator) fiuctuates around 
global mean {pt) = J2k=iPt/-^ev with probability distribution 

P{pi) = e-(^"*-<P*>)'/''^^-* (44) 

v2vr af. 



'Pt 

In this case Ylk=iiPt " {Pt)Y l^ev = o'lt and one obtains 

C2=-i= = c:^^^- + aj (45) 

Thus €2^'^'^ correlator contains (t|^ contribution from event-by-event fiuctuations of ob- 
servable mean pt and €2^'^ contribution from genuine two-particle correlations. This can be 
verified directly using a suitable MC event generator. For fiuctuations of observable mean 
given by Eq. (l44|) one has Yl,k=i^t ~ {Pt)Y/Nev ~ 3(t|j and in agreement with Eq. (l2T]) 

^^calc _ ^irne + 3 ^4 + g aJCa*"""^ (46) 

Assuming Cl^^^ one can try to separate cr|^ and C2 contributions as 

which is a solution of Eqs. (14511461) for vanishing 4-particle correlations {01^^^^= 0). 

IX. CONCLUSIONS 

A simple analytical calculation has shown that systematical shifts A„ in calculated values 
of correlators are generated if observable mean of the quantity under study is multiplicity- 
dependent x{nk) 7^ const. One can subtract such systematical effects from the calculated 
correlators C^^^^ to obtain "true" correlators using: C^™*' = C^^^^— A„. 

X. APPENDIX 

Fixed-multiplicity correction term given by Eq. fl32|) contains quantity (Jni^i'i^k) to be 
determined from MC simulation. We will show that width cr„(^,nfc) of n*°* fluctuations at 
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given fixed measured nl can be related to the width of fluctuations at fixed 



as (7ni^, TT-l) = ^ ■ c"nfc('^*°*) where ^ is acceptance parameter ^ = n*°Y (n^) and the width a, 
can be expressed as <Jn^{n]^^) = {(T)^/ri^. 

To demonstrate this let us divide given detector acceptance Vt^ into acceptance sub- 
regions oji + UJ2 + . . . + uj]^ = in such a way that in every acceptance region Ui equal 
average number of particles {n'^^) = {h^^) = {ni^)/N will be measured. Using N = fik 
acceptance subregions one has (ri'^') ^ 1 and Y^^=i{n^") = {nk)- For events with fixed total 
multiplicity n*°* and average measured multiplicity (n^) the number of particles n^' observed 
in acceptance region Ui will fluctuate (event-by-event) with some probability distribution 
Puji{n%^) characterized by the mean {h^^) and variance o^;. . 

Assuming that probability distributions P^^-{n^^) satisfy conditions for the applicability 
of generalized (e.g. m-dependent) Central Limit Theorem one can write 

Pn,{nk) = PnAT^^k ~ = ^n-(^fc) = P{nkK') (48) 



i=l 



where cr^, ^ +< + ••• + < = ^^^^ (denoting {Eli^^jN)'^' = (^))- This 

suggests that for (fik) = n*"*/^ large enough the probability distribution of events with 
measured multiplicities Uk in the set of events with fixed total multiplicity tends to be 
Gaussian and its width an,, increases with measured multiplicity as 



(Jn^=C-^ {Uk) (49) 

One can ask a similar question in the other way around: What are the fluctuations of 
total multiplicity of particles for events with fixed observed multiplicity fik ? 

Using the Bayes' theorem Gj one can calculate probability P(n*°*|njt) of observing the 
event with total multiplicity in the subset of events with fixed measured multiplicity hk 

1'''=^ - /P(n,|nt°t)P(nt°t)rfnt°t 
where P{A\B) denotes a conditional probability of observing quantity A for given B and 
P(ri*°*) is the prior probability^ of event with total multiplicity n*°*. Choosing the sample of 
events with P(?2*°*) = const simplifies the situation and one has P(n*°*|nfc) = A ■ P(nfc|n*°*) 
where A is a normalization constant. 



Note, that P(n*°') depends on a particular setting of the detector trigger. 
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For (fifc) large enough P(nfc|n*°*) given by Eq. (H8|) can be used to obtain probability 
P{n^°^\nk) of having event with total multiplicity in the group of events with fixed 
measured multiplicity -hk- For event sample with probability P(n*°*) = const one has: 

P{n'°'\nk) = X ■ P{fik\n'°') = X ^= = ^= (51) 

where we have denoted ^ ■ 0"^^^. = cxtot, used n*°Y^ = (n^) and normalized the resulting 
Gaussian distribution to unity. For event sample with probability P(n*°*) = 2u{l — un^°^) 
evaluation of denominator in Eq. (l50|) gives 2 ^ci}(l — ^ ri^ cj) and Eq. flSOj) then yields 

which is valid for fluctuations of ?t,*°* (at given fixed hk) within the range of P(n*°'^) distri- 
bution. For P{n"°^) = 2uj{l — cjn*"*) this means hkC, + 3crtot < 1/^^ which keeps ^hkCo < 1. 

For Poissonian probability P{hk) = {hk)^'' ■ e~^"''''^ /hk\ of measuring fik particles in a 
detector (exposed to events with fixed multiplicity n^°^=^{hk)) Bayes' law gives asymmetrical 
probability of n*°* for a given fixed hk'- 

r(n''"\hk) = X ^ . P(n'°') = X' ■ e-" ■ (n*°*)"^ ■ P(n'°') (53) 

Hkl 

(here A and A' are normalization factors). Fluctuations of at small fixed hk can be 
approximated by this function (see Fig.3 for V{n^°^\hk) at hk = 10 using P(n*°*) = const). 
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